A comparison of five multi attribute utility instruments

Abstract This paper presents the results of the validation study carried out to evaluate the Assessment of Quality of Life (AQoL) Instrument for the measurement of health related quality of life and utility. It involves, inter alia, the largest comparison of utility instruments that has been carried out to date. The five instruments included in the study are the AQoL, the Canadian HUI III, the Finnish 15D, the EuroQoL (EQ5D) and the SF36 with UK utility weights as quantified by Brazier (1998). The paper compares: (i) the absolute utility score obtained by different sub-populations; (ii) instrument sensitivity; (iii) the incremental differences in utility between different health states; (iv) the structural properties of descriptive systems; and (v) a limited comparison with a Time Trade-Off (TTO) assessment of own health by individuals. Using these criteria the AQoL performs very well. Its predicted utilities are very similar to those obtained from the HUI. There is evidence that the AQoL has greater sensitivity to health states than other instruments and its psychometric properties, as usually judged, are excellent. Despite this, it is concluded that, at present, no single MAU system can claim to be the gold standard and that researchers should select an instrument that is sensitive to the health states which they are investigating and that caution should be exercised in treating any of the instrument results as representing a utility score which truly represents a trade-off between life and health related quality of life

Static and lattice vibrational energy differences between polymorphs

A computational study of 1061 experimentally determined crystal structures of 508 polymorphic organic molecules has been performed with state-of-the-art lattice energy minimisation methods, using a hybrid method that combines density functional theory intramolecular energies with an anisotropic atom–atom intermolecular model. Rigid molecule lattice dynamical calculations have also been performed to estimate the vibrational contributions to lattice free energies. Distributions of the differences in lattice energy, free energy, zero point energy, entropy and heat capacity between polymorphs are presented. Polymorphic lattice energy differences are typically very small: over half of polymorph pairs are separated by less than 2 kJ mol?1 and lattice energy differences exceed 7.2 kJ mol?1 in only 5% of cases. Unsurprisingly, vibrational contributions to polymorph free energy differences at ambient conditions are dominated by entropy differences. The distribution of vibrational energy differences is narrower than lattice energy differences, rarely exceeding 2 kJ mol?1. However, these relatively small vibrational free energy contributions are large enough to cause a re-ranking of polymorph stability below, or at, room temperature in 9% of the polymorph pair

Modelling of crystal structure of cis-1,2,3,6 and 3,4,5,6-tetrahydrophthalic anhydrides using lattice energy calculations

Lattice energy calculations using a model potential were performed to model the crystal structures of cis-1,2,3,6- and 3,4,5,6-tetrahydrophthalic (THP) anhydrides. The optimized molecular models using the DFT method at the B3LYP/6-31G** level were found consistent with the available experimental evidence and allowed all differences observed in crystal packing between cis-1,2,3,6- and 3,4,5,6-THP anhydrides to be reproduced. Calculations provide evidence for the presence of dipole–dipole C=O?C=O intermolecular interactions and support the idea that the molecules distort from their ideal geometries, improving packing in both crystals. The search for minima in the lattice energy of both crystals amongst the more common space groups with Z’?=?1, using a simulated annealing crystal structure prediction procedure followed by lattice energy minimization showed that the observed structure of 3,4,5,6-THP anhydride (Z’?=?2) is the thermodynamically most stable, and allowed us to justify why 3,4,5,6-THP anhydride crystallizes in such a complex structure with 16 molecules in the unit cell. The computational model was successful in predicting the second observed form at 173 K for cis-1,2,3,6-THP anhydride as a polymorph, and could predict several hypothetical structures with Z’?=?1 that appear competitive with the observed structures. The results of phonon estimates of zero point intermolecular vibrational energy and entropy suggest that crystal structures of cis-1,2,3,6-THP anhydride cannot be predicted solely on the basis of lattice energy; factors other than thermodynamics favor the observed structures

Which conformations make stable crystal structures? Mapping crystalline molecular geometries to the conformational energy landscape

The ability to anticipate the shape adopted by flexible molecules in the solid state is crucial for engineering and predicting crystal packing and, hence, properties. In this study, the conformations adopted by flexible molecules in their crystal structures are assessed in terms of their relationship to the calculated global conformational landscape. The study quantifies the limits on molecular strain that can be induced by intermolecular interactions in single-component crystal structures of molecules with no intramolecular hydrogen bonding, demonstrating that some molecules are distorted by up to 20 kJ/mol by crystal packing forces. Furthermore, we find that crystallisation often selects high energy conformers, but only when the high energy conformer is more extended than the lower energy options, allowing for greater intermolecular stabilisation. Based on these observations, we propose that the crystallisability of conformers is assessed in terms of their energies and surface areas. We formulate this as a parameterised pseudo-energy related to molecular surface area, which leads to a dramatic improvement in our ability to predict the conformations adopted by molecules in their crystal structure

Accurate force fields and methods for modelling organic molecular crystals at finite temperatures

We present an assessment of the performance of several force fields for modelling intermolecular interactions in organic molecular crystals using the X23 benchmark set. The performance of the force fields is compared to several popular dispersion corrected density functional methods. In addition, we present our implementation of lattice vibrational free energy calculations in the quasi-harmonic approximation, using several methods to account for phonon dispersion. This allows us to also benchmark the force fields' reproduction of finite temperature crystal structures. The results demonstrate that anisotropic atom-atom multipole-based force fields can be as accurate as several popular DFT-D methods, but have errors 2-3 times larger than the current best DFT-D methods. The largest error in the examined force fields is a systematic underestimation of the (absolute) lattice energy

An optimized intermolecular force field for hydrogen bonded organic molecular crystals using atomic multipole electrostatics

This is the author accepted manuscript. It is currently under an indefinite embargo pending publication by the International Union of Crystallography.We present a re-parameterization of the a popular intermolecular force field for describing intermolecular interactions in the organic solid state. Specifically, we optimize the performance of the exp-6 force field when used in conjunction with atomic multipole electrostatics. We also parameterize force fields that are optimized for use with multipoles derived from polarized molecular electron densities, to account for induction effects in molecular crystals. Parameterization is performed against a set of 186 experimentally determined, low temperature crystal structures and 53 measured sublimation enthalpies of hydrogen bonding organic molecules. The resulting force fields are tested on a validation set of 129 crystal structures and show improved reproduction of the structures and lattice energies of a range of organic molecular crystals compared to the original force field with atomic partial charge electrostatics. Unit cell dimensions of the validation set are typically reproduced to within 3% with the re-parameterized force fields. Lattice energies, which were all included during parameterisation, are systematically underestimated when compared to measured sublimation enthalpies, with mean absolute errors of between 7.4 and 9.0%.Engineering and Physical Sciences Research Council (Doctoral Training Account studentship

An optimized intermolecular force field for hydrogen bonded organic molecular crystals using atomic multipole electrostatics

We present a re-parameterization of the a popular intermolecular force field for describing intermolecular interactions in the organic solid state. Specifically, we optimize the performance of the exp-6 force field when used in conjunction with atomic multipole electrostatics. We also parameterize force fields that are optimized for use with multipoles derived from polarized molecular electron densities, to account for induction effects in molecular crystals. Parameterization is performed against a set of 186 experimentally determined, low temperature crystal structures and 53 measured sublimation enthalpies of hydrogen bonding organic molecules. The resulting force fields are tested on a validation set of 129 crystal structures and show improved reproduction of the structures and lattice energies of a range of organic molecular crystals compared to the original force field with atomic partial charge electrostatics. Unit cell dimensions of the validation set are typically reproduced to within 3% with the re-parameterized force fields. Lattice energies, which were all included during parameterisation, are systematically underestimated when compared to measured sublimation enthalpies, with mean absolute errors of between 7.4 and 9.0%

Co-crystallisation of cytosine with 1,10-phenanthroline: computational screening and experimental realisation

Attempts to co-crystallise the nucleobases adenine, thymine, guanine, and cytosine with 1,10-phenanthroline by ball milling and solvent evaporation methods are described. A 1:1 co-crystal of cytosine and 1,10-phenanthroline can be obtained by grinding or by solvent evaporation. The structure contains two crystallographically independent cytosine and two independent 1,10-phenanthroline molecules (Z′ = 2). The cytosine molecules form two similar but crystallographically independent hydrogen-bonded chains, while the 1,10-phenanthroline molecules are arranged in π-stacks. Between the chains of cytosine and the π-stacks exist N-H⋯N and C-H⋯N interactions. Crystal structure prediction (CSP) calculations were applied to all four systems to assess their potential for co-crystallisation as well as the likely structures and intermolecular interactions that could result from co-crystallisation. Calculations on the cytosine system demonstrate that co-crystallisation results in a lower energy than the crystalline forms of the two starting materials, in line with the co-crystal formation observed. For the systems which did not form a co-crystal, CSP was used to explore potential packing arrangements, but found none which were lower in energy than that of the pure crystalline forms. In these cases there is significant disruption to the nucleobase hydrogen bonding between the pure compound and the hypothetical co-crystal. For pure adenine and guanine, the hydrogen-bonded ribbons form sheets which must be broken, whereas for thymine, the lack of hydrogen bond donors does not allow the hydrogen bonding present for pure thymine to be maintained while forming thymine-1,10-phenanthroline hydrogen bonds

Co-crystallisation of cytosine with 1,10-phenanthroline: computational screening and experimental realisation

Attempts to co-crystallise the nucleobases adenine, thymine, guanine, and cytosine with 1,10-phenanthroline by ball milling and solvent evaporation methods are described. A 1 : 1 co-crystal of cytosine and 1,10-phenanthroline can be obtained by grinding or by solvent evaporation. The structure contains two crystallographically independent cytosine and two independent 1,10-phenanthroline molecules (Z? = 2). The cytosine molecules form two similar but crystallographically independent hydrogen-bonded chains, while the 1,10-phenanthroline molecules are arranged in ?-stacks. Between the chains of cytosine and the ?-stacks exist N–H?N and C–H?N interactions. Crystal structure prediction (CSP) calculations were applied to all four systems to assess their potential for co-crystallisation as well as the likely structures and intermolecular interactions that could result from co-crystallisation. Calculations on the cytosine system demonstrate that co-crystallisation results in a lower energy than the crystalline forms of the two starting materials, in line with the co-crystal formation observed. For the systems which did not form a co-crystal, CSP was used to explore potential packing arrangements, but found none which were lower in energy than that of the pure crystalline forms. In these cases there is significant disruption to the nucleobase hydrogen bonding between the pure compound and the hypothetical co-crystal. For pure adenine and guanine, the hydrogen-bonded ribbons form sheets which must be broken, whereas for thymine, the lack of hydrogen bond donors does not allow the hydrogen bonding present for pure thymine to be maintained while forming thymine-1,10-phenanthroline hydrogen bonds

Construction of the descriptive system for the assessment of quality of life AQoL-6D utility instrument

BackgroundMulti attribute utility (MAU) instruments are used to include the health related quality of life (HRQoL) in economic evaluations of health programs. Comparative studies suggest different MAU instruments measure related but different constructs. The objective of this paper is to describe the methods employed to achieve content validity in the descriptive system of the Assessment of Quality of Life (AQoL)-6D, MAU instrument.MethodsThe AQoL program introduced the use of psychometric methods in the construction of health related MAU instruments. To develop the AQoL-6D we selected 112 items from previous research, focus groups and expert judgment and administered them to 316 members of the public and 302 hospital patients. The search for content validity across a broad spectrum of health states required both formative and reflective modelling. We employed Exploratory Factor Analysis and Structural Equation Modelling (SEM) to meet these dual requirements.Results and DiscussionThe resulting instrument employs 20 items in a multi-tier descriptive system. Latent dimension variables achieve sensitive descriptions of 6 dimensions which, in turn, combine to form a single latent QoL variable. Diagnostic statistics from the SEM analysis are exceptionally good and confirm the hypothesised structure of the model.ConclusionsThe AQoL-6D descriptive system has good psychometric properties. They imply that the instrument has achieved construct validity and provides a sensitive description of HRQoL. This means that it may be used with confidence for measuring health related quality of life and that it is a suitable basis for modelling utilities for inclusion in the economic evaluation of health programs.<br /